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Abstract. We consider the generalized evolution of compact level sets by 
functions of their normal vectors and second fundamental forms on a Rie- 
mannian manifold M. The level sets of a function u : M — > K evolve in such a 
way whenever u solves an equation Ut + F(Du, D 2 u) = 0, for some real func- 
tion F satisfying a geometric condition. We show existence and uniqueness of 
viscosity solutions to this equation under the assumptions that M has non- 
negative curvature, F is continuous off {Du = 0}, (degenerate) elliptic, and 
locally invariant by parallel translation. We then prove that this approach is 
geometrically consistent, hence it allows to define a generalized evolution of 
level sets by very general, singular functions of their curvatures. For instance, 
these assumptions on F are satisfied when F is given by the evolutions of 
level sets by their mean curvature (even in arbitrary codimension) or by their 
positive Gaussian curvature. We also prove that the generalized evolution is 
consistent with the classical motion by the corresponding function of the cur- 
vature, whenever the latter exists. When M is not of nonnegative curvature, 
the same results hold if one additionally requires that F is uniformly continu- 
ous with respect to D 2 u. Finally we give some counterexamples showing that 
several well known properties of the evolutions in R n are no longer true when 
M has negative sectional curvature. 



1. Introduction 

In the last 30 years there has been a lot of interest in the evolution of hypersur- 
faces of R™ by functions of their curvatures. In this kind of problem one is asked 
to find a one parameter family of orientable, compact hypersurfaces T t which are 
boundaries of open sets Ut and satisfy 

(1.1) V = -G(v, Dv) for t > 0, x e r t , and 

r t | t=0 = r 

for some initial set Tq = dUo, where v is the normal velocity of Tt, v = v(t, ■) is a 
normal field to T t at each x, and G is a given (nonlinear) function. 

Two of the most studied examples are the evolutions by mean curvature and 
by (positive) Gaussian curvature. In both cases, short time existence of classical 
solutions has been established. For strictly convex initial data Uq, it has been shown 
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that Ut shrinks to a point in finite time, and moreover, Tt becomes spherical at the 
end of the contraction. See [3l [TBI [2lJ 1221 Ell E3 [37] and the references therein. 

For dimension n > 3 it has been shown [20) that a hypersurface evolution T t may 
develop singularities before it disappears. Hence it is natural to try to develop weak 
notions of solutions to (jl.ip which allow to deal with singularities of the evolutions, 
and even with nonsmooth initial data Tq. 

There are two mainstream approaches concerning weak solutions of (|1.1|) : the 
first one uses geometric measure theory to construct (generally nonuniquc) varifold 
solutions, see [51(21], while the second one adapts the theory of second order vis- 
cosity solutions developed in the 1980's (see [9] and the references therein) to show 
existence and uniqueness of level-set weak solutions to ([l.ip . 

In this paper we will focus on this second approach. The first works to develop 
a notion of viscosity level set solution to (|1.1|) were those of Evans and Spruck 
[13] and, independently developed, Chen, Giga and Goto [7J. This was followed by 
many important developments, which we find impossible to properly quote here; 
we refer the reader to the very comprehensive monograph [17j and the bibliography 
therein. This level set approach consists in observing that a smooth function u : 
[0, T] x R" — > R with Du := D x u ^ has the property that all its level sets evolve 
by ([l.ip if and only if u is a solution of 

(1.2) u t + F(Du,D 2 u) = 0, 

where F is related to G in (|1.1|) through of the following formula: 

d-3) FM = ^ G {k¥\{ I - P -w) A 

The function F is assumed to be continuous off {p = 0} and (degenerate) elliptic, 
that is 

(1.4) F(p,B) < F(p,A) whenever A < B. 
Because of (|1.3p . F also has the following geometric property: 

(1.5) F(Xp, XA + up <g> p) = XF{p,A) for all A > 0,/zG R. 

2 

The function F does not generally admit any continuous extension to R™ x R™ 
but, if it is bounded near {p — 0} (this is the case of the mean curvature evolution 
equation), one can show that there is a unique viscosity solution to (|1.2p with initial 
datum u(0, x) = g(x) (for any continuous g such that Tq — {x : g(x) = 0}). Next 
one can also see that if 9 : R — ► K. is continuous and u is a solution of (|1.2p then 
9 o u is a solution too, and this, together with a comparison principle, allows to 
show that the generalized geometric evolution 

T -> r t := {x : u(t, x) = 0} 

is well defined (that is, the zero level set of a solution to (|1.2p only depends on the 
zero level set of its initial datum). It is also possible to show that this level set 
evolution agrees with any classical solution of (jl.ip . 

When F(p, A) is not bounded as p — > (this is the case of more singular equations 
such as the Gaussian evolution), then the standard notion of viscosity solution to 
(II. 2p (as is used, for instance, in [7J) at points z = (t,x) where the test function ip 
satisfies Dip(z) = is not suitable to tackle the problem. In this case two different 
modifications of the notion of solution have been proposed in the literature. 
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One possibility is simply not to specify any condition for the derivatives of a test 
function tp such that u — tp attains a maximum or a minimum at a point (to,xo) 
with D<p(to,xo) — 0- This is Goto's approach in [T^]. When one uses this definition 
of solution, the corresponding comparison theorem becomes harder to prove, and 
it is indeed a stronger statement since the class of solutions becomes bigger in this 
case, while the existence result is comparatively weaker. 

The other possibility is to make the class of test functions tp smaller, in a clever 
way so that, if Zk — > zq and Dtpk(zk) — * 0, one can show that F(Dtpk(zk), D 2 tpk{zk)) 
goes to 0, and then to demand that a subsolution u should satisfy that if u— tp has a 
maximum at zq then <ft(zo) < 0. This is what Ishii and Souganidis did in |30j . The 
corresponding (sub)solutions are called T- (sub) solutions in Giga's book [17] . In 
this approach the maximum principle is relatively easier to prove, while existence 
becomes harder (and is really a stronger result, because the class of solutions is 
smaller in this case). 

The aim of this paper is to investigate to what extent one can develop a general 
theory of (viscosity) level-set solutions to the problem of the evolution of hypersur- 
faces by functions of their curvatures in a Riemannian manifold. To the best of our 
knowledge, the only work in this direction is Ilmanen's paper [26] (in fact this is 
the only paper we know of in which second order viscosity solutions are employed 
to deal with a second order evolution equation within the context of Riemannian 
manifolds). In [26 Ilmanen shows existence and uniqueness of a (standard) vis- 
cosity solution to the mean curvature evolution equation, that is (JTTTJ) in the case 
when F is given by 

FfeA) = -trace ((j-^)^), 

with initial condition u(0, x) — g(x), thus obtaining a corresponding generalized 
evolution by mean curvature, some of whose geometric properties he next studies. 
For instance, he proves that if noncompact initial data T are allowed then one 
loses uniqueness of the generalized geometric evolution. 

In recent years, an interest has grown in the use of viscosity solutions of (first 
order) Hamilton- Jacobi equations defined on Riemannian manifolds (in relation to 
dynamical systems, to geometric problems, or from a theoretical point of view), see 
[32l HU 021 HI [311 HOI E3] , but no second order theory, apart from Ilmanen's paper, 
has apparently been developed for parabolic equations (in the case of stationary, 
degenerate elliptic equations, such a study was recently started in [5]). 

We believe that a level set method for generalized evolution of hypersurfaces by 
functions of their curvatures can be useful in the setting of Riemannian manifolds. 
On the one hand we think that it is very natural, from a geometric point of view, 
to try to study the evolutions of level sets in a general Riemannian manifold M 
by their Gaussian (or by other functions of their) curvatures, in a way that is 
supple enough so that nonsmooth initial data and singularities of the evolutions 
are allowed. On the other hand, as one sees, for instance, by restricting to the case 
M = K" endowed with a non Euclidean metric, the tools developed here allow to 
treat level set evolutions in inhomogeneous media, in which the function F depends 
(in a very special manner) on the position variable x. 
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Let us briefly describe the main results of this paper. In Section 2 we consider 
equations of the form (jl.ip . (|1.2|) for level sets of functions u defined on a Rie- 
mannian manifold M, and we show how the F's corresponding to the evolutions by 
mean curvature (even in arbitrary codimension, in the line of [2]) and by (positive) 
Gaussian curvature are extended to Jq(M) in such a way that F is (degenerate) 
elliptic, translation invariant, geometric, and continuous off {Du = 0} (see prop- 
erties (A - D) in Section 2 below). Following j30j[T7], for each F we next define 
an appropriate class of test functions A{F) which allows us to deal with equation 
(II. 2[) on M, and we define the corresponding class of ^"-solutions, see Definitions 
12. 4[ 12.71 below. We also show that for all F which are continuous off {Du = 0}, 
elliptic, translation invariant and geometric, one has that A(F) ^ provided that 
M is compact. Moreover, in the cases when F is given by the mean curvature or 
the Gaussian curvature evolution equations, we have A{F) ^ no matter whether 
M is compact or not. 

In Section 3 we present some technical results that will be used later on in the 
proofs of the main results. 

Section 4 is devoted to proving a comparison result for viscosity solutions of (|1.2p 
on M: under the above assumptions on F (namely, continuity, ellipticity, geometric- 
ity and translation invariance) we show that if M has nonnegative curvature u is a 
subsolution, v is a supersolution, u < v on {0} x M, and limsup( t x )_ (00 (w — v) < 
(this condition is understood to be requiring nothing when M is compact), then 
u < v on [0, T] x M. When M is not of nonnegative curvature, we have to addi- 
tionally require that F be uniformly continuous with respect to D 2 u. 

In Section 5 we show that Perron's method (first used in [28 ) works to produce 
^4(-F)-solutions of (|1.2p on a Ricmannian manifold M, provided that comparison 
holds and A(F) ^ 0. 

Therefore, for all such M and F, for every compact subset Tq of M, and for 
every continuous function g on M such that Tq = {x M : g(x) = 0}, there exists 
a unique solution of (|1.2p on M with initial condition u(0, •) = g. One can then 
define, for each compact Tq, an evolution T t = {x G M : u(t,x) — 0}, t > 0. In 
Section 6 we see that T t does not depend on the function g chosen to represent Tq, 
and consequently the generalized geometric evolution Tq t— > T t is well defined. 

Next, in Section 7 we prove that this generalized evolution is consistent with 
the classical motion, whenever the latter exists. Namely, if (T t ) t ^ T , is a family of 
smooth, compact, orientable hypersurfaces in a Riemannian manifold M evolving 
according to a classical geometric motion, locally depending only on its normal 
vector fields and second fundamental forms according to an equation of the form 
(jl.ip . and Tq can be represented as the zero level set of a smooth function g on M, 
then T t coincides with the generalized level set evolution (with initial datum Tq) 
defined above. 

Finally, in Section 8, we give some counterexamples showing that several well 
known properties of generalized solutions to the mean curvature flow cannot be 
extended from Euclidean spaces to Riemannian manifolds of negative sectional 
curvature. For instance, Ambrosio and Soner [2] showed that the distance function 
from T t C M™ given by \d\(t, x) — dist(x,r t ) is a supersolution of (|1.2p when F 
corresponds to the mean curvature evolution equation. We show that this result fails 
when R ra is replaced with a manifold of negative sectional curvature. On the other 
hand, if M has negative curvature, then equation (|1.2p does not preserve Lipschitz 
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properties of the initial data, in contrast with [17\ Chapter 3] . And, again in the 
case of the mean curvature flow, if To, Tq are smooth 1-codimensional submanifolds 
of a manifold M of negative curvature, then the function t t— > dist(r t ,r t ) can be 
decreasing, in contrast with [13[ Theorem 7.3]. 

An the end of this article, the reader will find an appendix describing a compar- 
ison and an existence result for (standard) viscosity solutions to general evolution 
equations of the form 

u t + F(x, t, u, Du, D 2 u) = 

where F has no singularities. We omit the proofs because they resemble (and are 
easier than) those of the main comparison and existence result for ^-solutions of 
(|1.2|) given in Sections 4 and 5. 

Notation. M will always be a finite-dimensional Riemannian manifold. We will 
write (■, ■) for the Riemannian metric and | ■ | for the Riemannian norm on M. 
The tangent and cotangent space of M at a point x will be respectively denoted 
by TM X and TM*. We will often identify them via the isomorphism induced 
by the Riemannian metric. The space of bilinear forms on TM X (respectively 
symmetric bilinear forms) will be denoted by £ 2 (TM X ) (resp. C 2 (TM X )). Elements 
of £ 2 {TM X ) will be denoted by the letters A,B,P,Q, and those of TM* by 
etc. Also, we will respectively denote the cotangent bundle and the tensor bundle 
of symmetric bilinear forms in M by 

TM* := (J TM*, T 2 , S (M) := (J C 2 S (TM X ). 

xEM x£M 

We will also consider the two-jet bundles: 

J 2 M := |J TM* x C 2 S (TM X ), J 2 (M) := |J (TM* \ {0 X }) x C 2 S (TM X ). 

The letters A, Y, Z will stand for smooth vector fields on M, and VyX will always 
denote the covariant derivative of A along Y. Curves and geodesies in M will be 
denoted by 7, c, and their velocity fields by j',a'. If A is a vector field along 7 
we will often denote X'(t) = ^X(t) — Vy( t )A(i). Recall that A is said to be 
parallel along 7 if X'(t) = for all t. The Riemannian distance in M will always 
be denoted by d(x, y) (defined as the infimum of the lengths of all curves joining x 
to y in M). 

Given a smooth function u : M — > R, we will denote its differential by D x u £ 
TM*; its gradient vector field will be written as Vit, and its Hessian as D 2 u. Recall 
that, for any two vector fields A, Y satisfying A(p) = v, Y(p) = w at some p € M 
we have: 

D 2 x u (A, Y) := (V F V«, A) , D 2 x u (v, w) := D 2 x u (A, Y) (p). 
Given a function v : M — > M we will use the notation: 

v* (t, x) = lim ri0 sup{v(s, y) : y <G M, s > 0, \t - s\ < r, d(y, x) < r}, 
v*(i, x) — lim r | inf{u(s, y) : y G M, s > 0, \t — s\ < r, d(y 7 x) < r}; 

that is v* denotes the upper semicontinuous envelope of v (the smallest upper 
semicontinuous function, with values in [—00,00], satisfying v < v*), and similarly 
v* stands for the lower semicontinuous envelope of v. 

We will make frequent use of the exponential mapping exp^, and of the parallel 
translation along a geodesic 7. Recall that for every x € M there exists a mapping 
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exp x , defined on a neighborhood of in the tangent space TM X , and taking values in 
M, which is a local diffeomorphism and maps straight line segments passing through 
onto geodesic segments in M passing through x. The exponential mapping induces 
a local diffeomorphism on the cotangent space TM* , via the identification given by 
the metric, that will be also denoted by exp^,. On the other hand, for a minimizing 
geodesic 7 : [0, i] — > M connecting x to y in M, and for a vector v £ TM X there is a 
unique parallel vector field P along 7 such that P(0) = v, this is called the parallel 
translation of v along 7. The mapping TM X 3 v 1— > P(£) £ ^Af^ is a linear isometry 
from TjWj onto TM y which we will denote by L xy . This isometry naturally induces 
an isometry between the space of bilinear forms on TM X and the space of bilinear 
forms on TM y . Whenever we use the notation L xy we assume implicitly that x and 
y are close enough to each other so that this makes sense. 

By im(x) we will denote the injectivity radius of M at x, that is the supremum 
of the radius r of all balls B(0 X , r) in TM X for which exp^. is a diffeomorphism from 
B(Q x ,r) onto B(x,r). Similarly, i[M) will denote the global injectivity radius of 
M, that is i(M) = inf{ijw(:r) : x S M}. Recall that the function x 1— » im{x) is 
continuous. In particular, if M is compact, we always have i(M ) > 0. 

2. General curvature evolution equations on Riemannian manifolds 
Consider the following evolution equation on a Riemannian manifold M, given 

by 

(CEE) u t - F(Du,D 2 u) =0 on (0,T)xM, 

u(0,x) — g(x), on x S M, 

where u is a function of (t,x) £ [0, T) x M. 

In what follows, «t, Du and D 2 it will stand for D t u, D x u(t,x) £ TM* and 
D^u(t,x) £ C 2 S (TM X ), respectively. The function F is assumed to be continuous 
on the normal vector to the level set T t — {x £ M : u(t, x) = 0} and on the 
curvature tensor, and having the form 

(2.1) F( ^ ) = |C|G (^(,-^), 

for all C £ TM*\{0 x } and A £ £ 2 (TM X ), where G is any (nonlinear) function such 
that: 

(A) F : Jq (M) -> E is continuous; 

(B) F is (degenerate) elliptic, that is, 

A < B =>■ F({,B)<F(£,A), 

for all x £ M, £ G TM* \ {0}, A, B £ C 2 S {TM X ); 

(C) F is translation invariant, meaning that there exists r > such that 

F(L xy C,A)=F(C,L yx (A)), 
for every i, y£ M, d(x,y)<T, £ G TM* \ {0}, ^4 G C 2 (TM y ). 
Notice that, because 

(/-^) K«0 = o. 

any function F of the form ()2.ip also satisfies 
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(D) F is geometric, that is, 

F(\(,\A + fiC®() = \F(t,A) 
for every A > 0, fi € R. 

Two very important problems where such functions F arise are the evolutions 
of level sets by mean curvature and by Gaussian curvature. 

Example 2.1. Motion of level sets by their mean curvature. 

If u is a function on [0, T]x M such that Du(t, x) ^ for all t, x with u(t, x) — c, 
then each level set T t = {u(t, •) = c} evolves according to its mean curvature if and 
only if u satisfies 

u t ,. ( Du 
div 



|£>u| " \\Du\ 

(that is, the normal velocity of T t at a point x equals (n — 1) times the mean 
curvature of T t at x), which in turn is equivalent to 

Du ® Du^ 
~\Du~\ 



(MCE) u t - trace ( ( I , 2 ) D 2 u ) = on (0, T) x M. 



That is, u t + F(Du, D 2 u) = 0, where 
(2.2) F ( C7 ^) = _traccff/-^)A 



ICI 2 

It is not difficult to see that the function F : Jq(M) — ► R is continuous (though 
the function F remains undefined at £ = and, in fact, there is no continuous 
extension of F to J 2 (M). Nevertheless, F((,A) remains bounded as £ — » 0). 

Let us now check that the function F is degenerate elliptic. If P < Q, since 

c® c 

R:= I — pr- > and S := Q — P > 0, we obtain from the properties of the trace 

\c\ 

that trace(i?5) > and therefore 

(2.3) F(C, P) - F(C, Q) = trace ( (i - ^f) (Q - P) \ > 0. 



Finally, let us see that the function F in (|2.2[) is translation invariant. Notice 
that trace(A) = trace(L xy oio L xy ) = tiace(L yx (A)), and 

trace ( ^jtS ° Ly X (A) J = trace ( L xy o -r^»~ o L yx (A) o L^ 1 



|^|2 -~V-V"7 I ~~ \ - |£|2 ■"» ar v"7 - ^a:j/ 

C®C _ r _i 



- traCe l-^xi/ - |^| 2 " ~a:v 

On the other hand, 

{2A) Lxy ICI 2 l«l 2 ' 

hence we immediately deduce that F(L xy (, A) — F((, L yx (A)) whenever d(x,y) < 
i(M), C e TMj;, A e C 2 s (TM y ). 



Example 2.2. Motion of level sets by their Gaussian curvature. 
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Now, if u is a function on [0, T] x M such that Du(t, x) ^ for all t, x with 
u(t,x) = c, then all level sets T t — {u(t, •) = c} evolve according to their Gaussian 
curvature if and only if u satisfies 

Ut dct(W Vu 



\Du\ \ \\Vu\ 

where V T stands for the orthogonal projection onto TT t of the covariant derivative 
in M. This equation is equivalent to 

1 ( T Du®Du\ 2 Du®Du 



(GCE) ut ~ '^ |det 1 W I 7 - ~WT J 79 u + "W 1 = 1 

That is, u t + H(Du, D 2 u) = 0, where 

tf(C,A) = H C|dc«((,-^).4 + ^ 

However, the function H is not elliptic, so this problem cannot be treated, in its 
most general form, with the theory of viscosity solutions. Nevertheless, if our initial 
data u(0, x) = g(x) satisfies that D 2 g(x) > (that is, if the initial hypersurface T = 
{x € M : g(x) = c} has nonnegative Gaussian curvature) then it is reasonable, and 
consistent with the classical motion of convex surfaces by their Gaussian curvature, 
to assume that D 2 u(t,x) > for all (t,x) with u(t,x) = c (that is, r t will have 
nonnegative Gaussian curvature as long as it exists). In this case our equation 
becomes 

( + GCE) u t - |^.|det + ( m (I - j D 2 u + J = 0, 

where det + is defined by 

n 

det+(A) = TT max{Aj, 0} 

if Ai, A„ are the eigenvalues of A. That is, 

u t + F(Du,D 2 u) = 0, 

where 

(2.5) F(C ,A) = - K |det + af/-^V+ C0C 



.ICI V ICI 2 ) ICI 2 

As in the case of the mean curvature, it is not difficult to see that F is elliptic 
and translation invariant, and that F is continuous off {C = 0} (this time the 
singularities at £ = are of higher order, as F((, A) generally tends to ±oo as C 
goes to 0). 

Example 2.3. Motion by mean curvature in arbitrary codimension. 

If T is a fc-codimensional surface of an n-dimensional Riemannian manifold M, 
we choose a continuous function vq with To = u ( 7 1 (0), and consider 

u t + F(Du,D 2 u) = 0, u(0, x) = v (x), 

where 

d-k 



F((,A)=J2*i(Q) 
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and 

Ai(Q) < A 2 (Q) < ... < Xd-xiQ) 
are the eigenvalues of Q ;= S^AS^, with 



ICI 2 



corresponding to eigenvectors orthogonal to £ (note that £ is an eigenvector corre- 
sponding to the eigenvalue of Q). 

The same proof as in [2] shows that F is elliptic, the key observation is that 

\JQ) = maximin ^'J 1 ' : E C TM X , codim(S) < i - 1}. 
rieE \r]Y 

On the other hand, it is easy to see, as in Example [^T] above, that F is translation 
invariant. 



Our aim is to establish comparison, existence and uniqueness of viscosity solu- 
tions to the general curvature evolution equation ICEE1 and then to prove that the 
resulting generalized motion is consistent with the corresponding classical motion 
(whenever the latter exists). However, because this equation is, in general, highly 
singular, one has to define very carefully what a viscosity solution to ICEEI is at 
points where Du = 0. Here we will adapt Ishii-Souganidis' definition [3D] (see also 
[17) ) from the Euclidean to the Ricmannian setting. This requires a slight change 
in the definition of the set of test functions tp. 

Definition 2.4. Let F : Jg(M) — > R be continuous, (degenerate) elliptic, trans- 
lation invariant and geometric. Denote by T = T(F) the set of functions / £ 
C 2 ([0,oo)) such that /(0) = /'(0) = /"(0) = and f"{s) > for s > which 
satisfy 

(2.6) lim ^^(C, 21) = lim ^^F((, -21) = 0. 

ICHo |C| ICHo |C| 

It is clear that J 7 is a cone (that is, / + g 6 T and A/ € J- whenever f,g£j-,\£ 
[0,oo)). 

Proposition 2.5. If M is compact and F : Jq(M) — > R is continuous, elliptic, 
translation invariant, and geometric, then ^F(F) ^ 0. 

Proof. One can adapt the proof given in [3pJ p. 229] for the case M = W 1 . The only 
difference (apart from the replacement of / with 21) is that |C| = |CU depends on 
the point x such that £ £E TM X , and one has to be cautious about this dependence 
(as a matter of fact, that is why we require compactness of Al). Let us give the 
essential details for the reader's convenience. 

Since F is continuous on Jq(M) and the sets {(^,±2/) : Id,^ = l,x € M} are 
compact in Jq (M), there exists a continuous function c : (0, oo) — > (0, oo) such that 

-c(|d)<F(d2/)<F(C,-2/)< C (|d) 

for all C £ TM* \ {0 X : x £ M}. Without loss of generality one can then assume 
that c is C 1 on (0,oo) and satisfies (1/c)' > in (0,1], lim r ^ + c i r ) = °°> an d 
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lim r _^o+ (V c )'( r ) = 0- Then it is not difficult to show that an appropriate extension 
to [0, oo ) of the function / defined on [0, 1] by 

f(r) = { £ £) ds > if < r < 1; 
' \ 0, if r = 0, 

belongs to T(F). □ 

For many interesting choices of the function F it is easy to show that T{F) ^ 
without requiring M to be compact: 

Example 2.6. If F is given by (|2.2p (corresponding to the mean curvature evolu- 
tion equation), then we may take / € ^F{F) of the form 

f(t) = t\ 

On the other hand, when F is associated to the Gaussian curvature evolution 
equation (that is, F is given by Q2.5p ) then 

/(f) = t 2n 

belongs to T{F) (here n is the dimension of M). 

Definition 2.7. We define the set A(F) of admissible test functions for the equation 
(|CEE|) as the set of all functions ip € C 2 ((0,T) x M) such that, for every z = 
(t , x ) e (0, T) x M with Dip(z ) = there exist some S > 0, / € T, w e C([0, oo)) 
satisfying lim r .^ + w(r)/r = and 

\ip(z) - cp(z ) - tp t (z )(t - f )| < /(d(ar, ar )) + w(|< - f |) 

for all z = (t,x) € i3(zo7^)- 

Notice that in particular, for all <p e ^4.(-F) we have that 

£><p(z) = D 2 ip{z) = 0. 

Proposition 2.8. If M is a compact Riemannian manifold then the class A(F) of 
admissible test functions is dense in the space C(M) of continuous functions on M. 

Proof. It is not difficult to check that the class A(F) satisfies the hypotheses of the 
Stone- Weierstrass theorem. □ 

Definition 2.9. We will say that an upper semicontinuous function u : [0, T) x 
M — > R is a viscosity subsolution of (|CEEp provided that, for every ip S A(F) and 
every maximum point z = (t, x) of u — <p, we have 

U t + F(Dip(z),D 2 <p(z)) <0 ^ D(p(z)^0, 
\<Pt(z)<0 otherwise. 

Similarly, we will say that a lower semicontinuous function u : [0,T)xM^l is 
a viscosity supersolution of (|CEEp if, for every ip G A(F) and every minimum point 

z = (t, x) of u — <p, we have 

{p t + F(Dp{z),D 2 p(z)) >0 if^(z)^0, 
1 <Pt{ z ) > otherwise. 

A viscosity solution of (|CEEJ) is a continuous function u : [0, T) X M — ► K which is 
both a viscosity subsolution and a viscosity supersolution of (|CEEp . 
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In [17] this kind of solution is called an ^"-solution, but here we will simply call 
it a solution. It is clear that one can always assume that the minimum or maximum 
in these definitions are strict. 

Notice that the set of test functions ip we are using is smaller than the standard 
one in the general theory of viscosity solutions, and that we here require that ip is C 2 
with respect to the variables t and x (while in the usual definition of the parabolic 
semijets one demands C 1 differentiability with respect to t and C 2 differentiability 
with respect to x). 

It is easy to check that this definition is consistent with u being a classical 
solution. Indeed, if u is a classical solution then we have Du{z) ^ and Ut(z) + 
F(Du(z), D 2 u(z)) — for all z. Then, if <p G A(F) is such that u — <p attains a 
minimum at z, we have <pt(z) = Ut(z), Dip(z) = Du{z) ^ 0, and D 2 u(z) > D 2 ip{z). 
Since F is elliptic we get 

<p t (z) + F(Dp(z),D 2 (p(z)) > u t (z) + F(Du(z),D 2 u(z)) = 0, 

that is, u is a supersolution at z. A similar argument shows that u is a subsolution. 

It can be proved, as in the Euclidean case [17] . that if the lower and upper 
semicontinuous envelopes of F (denoted by F_ and F respectively) are finite and 
£(0,0) = F(Q, 0) = 0, then every standard viscosity solution is an ^"-solution, and 
conversely. This is the case of the F associated to the mean curvature evolution. 



3. Some technical tools 

In this section we collect some rather technical results that will be needed in the 
proof of the main comparison theorem. 

First, we will need to use the following variant of the maximum principle for 
semicontinuous functions already used in [5], which we restate here for the reader's 
convenience. 

Theorem 3.1. Let M%, ...,M% be Riemannian manifolds, and f2j C Mj open sub- 
sets. Define O = fli x . . .xfi„ C M\ x . . .xMk = M . Let Ui be upper semicontinuous 
functions on Qi, i = 1, k; let ip be a C 2 smooth function on ft and set 

u>(x) = Wi(xi) + . . . + u n (x k ) 

for x — (ari, ...,£&) G f2. Assume that (x%, . . . ,x k ) is a local maximum of lu — (p. 
Then, for each e > there exist bilinear forms Bi G ^((TMi)^), i = 1, k, such 
that 

2 -j- 

(D Xi (p(x),Bi) G J ' u.i(xi) 
for i = 1, k, and the block diagonal matrix with entries Bi satisfies 

f B, ... 
- + \\A\\\l< : •.. : (j A + cA-. 
V ... B k 
where A = D 2 <p(x) G £ 2 S (TM & ). 
We recall that 

J 2 + f(x) = {(Dtp(x), D 2 ip(x)) : ip G C 2 (M,R), f—ip attains a local maximum at x}, 
and 

J 2 - + f(x) = {(0 A) G TM* X x C S {TM X ) : 3(x k ,( k ,A k ) G M x TM* x C s {TM Xk ) 



12 D. AZAGRA, M. JIMENEZ-SEVILLA, F. MACIA 

s.t. (( k ,A k )eJ 2 -+f(x k ), (x k ,f(xk)Xk,A k )^ (x,f(x)X,A)}, 

see [5]. 

Another important ingredient of the proof of our main comparison result is the 
following Proposition, established in 5, Proposition 3.3]. 

Proposition 3.2. Consider the function ^(x,y) = d(x,y) 2 , defined on M x M. 
Assume that the sectional curvature K of M is bounded below, say K > —Ko. Then 

Dl^{x,y){v,L xv vf < 2K Q d(x iy ) 2 \\v\\ 2 

for all v £ TM X and x, y £ M with d(x, y) < mia{iM{x), *m(2/)}- 

In particular, if —Ko > (that is M has nonnegative sectional curvature) one 
has that the restriction of D 2 y ^/(x,y) to the subspace T> — {(y,L xy v) : v £ TM X } 
of TM X x TM y is negative semidefinite. 

We will also need the following auxiliary result. 

Lemma 3.3. Let 4> £ USC(M), ip £ LSC(M), f £ T{F), and 

m a := sup (<p(x) - ip(y) - af (d(x,y) 2 )) 

M x M 

for a > 0. Suppose m a < oo for large a and let {x a , y a ) be such that 
lim (m a - ((f>(x a ) - ip(y a ) - af{d(x a ,y a ) 2 ))) = 0. 

Then we have: 

(1) lim^oo af{d(x a ,y a ) 2 ) =0, and 

(2) if x £ M is a limit point of x a as a — > oo then 

lim m a = 4>{x) — ip(x) — sup (<fi(x) — ip(x)). 

Proof. A more general form of this result is proved in .9, Theorem 3.7] in the case 
when M is an Euclidean space, and the same proof clearly works in a general metric 
space. □ 

Let us now define V 2 ' + , V 2 ~ , V ' , and V ' , the "parabolic" variants of the 

2 -4- 2 

semijets J 2,+ , J 2 '~ , J , J introduced in [5] for functions defined on a Rie- 
mannian manifold. 

Definition 3.4. Let / : (0, T) x M — > (— oo, +oo] be a lower semicontinuous (LSC) 
function. We define the parabolic second order subjet of / at a point (to,xo) £ 
(0, T) x M by 

"P 2 '~f(to,%a) ■={(D t cp(to,Xo),D x cp(to,Xo)),Dlip(to,x )) : <p is once continuously 

differentiable in t £ (0,T), twice continuously differentiable in x £ M 
and / — if attains a local minimum at (to,xo)}- 

Similarly, for an upper semicontinuous (USC) function / : (0,T)xM-> [— oo, +oo), 
we define the parabolic second order superjet of / at (to,Xo) by 

^ 2 ' + /(io,^o) ■={(Dtip(to,Xo),D x ip(to,Xo)),Dlip(to,x )) : <p is once continuously 

differentiable in t £ (0, T), twice continuously differentiable in x £ M 
and / — tp attains a local maximum at (to,xo)}. 
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Observe that V 2 '-f(t,x) and V 2 >+f(t,x) are subsets of 1 x TM* x C 2 S (TM X ). 
Notice that we can assume that the auxiliary functions ip arc defined on a neigh- 
borhood of (to,Xo). We may as well assume (just by adding a function of the form 
±sd(x, xo) 4 ) that the minima or maxima in these definitions are strict. It is also 
easily seeing that the min or max can always be supposed to be global. 

Definition 3.5. Let / : (0, T) x M — > (— oo,+oo] be a LSC function and (t,x) € 
(0, T) x M, We define V 2 ^ f{t, x) as the set of the (a, (,i)eRx TM* x C 2 JTM X ) 
such that there exist a sequence (xk, Ofc, Cfci Ak) in M x E x TM* k x C 2 {TM Xh ) 
satisfying: 

i) K,a,Afc) eV 2 '-f{t k ,x k ), 

ii) limfe(4,a;fe,/(4,^fe),afe,Cfe,^fc) = {t,x, f(t,x),a, £,A). 

2 + 

The corresponding definition of V ' f(t,x) when / is an upper semicontinuous 
function is then clear. 

The next two lemmas are needed to establish the parabolic version of the maxi- 
mum principle we state below. 

Lemma 3.6 (|5J). Let U C M be an open subset, (t, z) £ (0, T) x U and a func- 
tion if : (0, T) x U — > K. Assume that if is once continuously differ •entiable in 
(0,T) and twice continuously differentiable in U. Define tJj(s,v) — ip(s.exp z v) on 
a neighborhood o/O £ TM Z . Let V be a vector field defined on a neighbourhood of 
in TM Z , and consider the vector field defined by V(y) = D exp z (w y )(V(w y )) on 
a neighbourhood of z, where w v := exp~ 1 (y), and let 

o y {r) = cxp z (w y + rV(w y )). 

Then we have that 

D 2 v i,{V, V)(t, w v ) = D 2 x f(V, V)(t, y) + (V^(i, y), <(0)). 
Observe that o~"(0) = so, when y = z, we obtain 

Dfy(t,0) = Dlip(t,z). 
Proof. Analogous to [5l Lemma 2.7]. □ 

Lemma 3.7. Let U C M be an open subset, (t, z) £ (0, T)xU and u : (0, T) x M —* 

[—oo, oo) be an upper semicontinuous function and consider a neighbourhood V of 
G TM Z and u : (0, T) x V — > [—00,00) defined as u(s,v) — u(s,exp z v). Then, if 
(b, C, A) e K x TM* x C 2 {TM Z ), 

(b,(,A)EV 2 ' + u(t,z) ^ (fc,C,i)eP 2 ' + i(t,0). 
Proof. Use the above Lemma as in the proof of Proposition 2.8]. □ 

As in [26] in the case of the mean curvature evolution equation and [5] in the 
case of general (nonsingular) stationary equations, the following result is one of 
the keys to the proof of the comparison result for general (nonsingular) evolution 
equations which we give in the Appendix. 
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Theorem 3.8. Let Mi, ...,Mk be Riemannian manifolds, and f2i C Mi open sub- 
sets. Define f2 = (0, T) X f2i X • • • X Qk- Let Ui be upper semicontinuous functions 
on (0,T) x Qi, i = 1, k; let ip be a function defined on Q such that it is once 
continuously differ entiable in t £ (0,T) and twice continuously differ entiable in 
x := (a;i j Xk) £ Oi x ■ ■ ■ x 0^ and set 

tu(t,xi, ...,x k ) = Ui(t,xi) H h u k (t,x k ) 

for (t, xi, .... Xk) £ fi. Assume that (i, X\, . . . , Xk) is a maximum of u) — tp in Q. 
Assume, moreover, that there is an r > such that for every M > there is C > 
such that for i = 1, k, 

ai < C whenever (a i: Q, A) € Pm, Ui(t, Zj) 

d(xi,Xi) + \t-t\<T and \ Ui (t, Xi )\ + \&\ + \\Ai\\ < M. 

Then, for each e > there exist real numbers bi and bilinear forms Bi € 
C 2 s ((TMi) i = 1, fc, smc/i i/iai 

(bi,D Xi ip(t, xi,... ,x k ),Bi) € V 2 ^ Ui(t,Xi) 
for i = 1, ...,k, and the block diagonal matrix with entries Bi satisfies 

... 




< A + eA , 



Bk 



where A = Dltp(t,xi, ...,Xk) and b± + ■ ■ ■ + bk — ^-(t, x\, x k ). 

Proof. The result is proved in [9] for M, = M. ni , i = l,...,k. As in the sta- 
tionary case [5], we can reduce the problem to this situation by an adecuate 
composition with the exponential mappings. Let us give some details for com- 
pleteness. We may assume (by taking smaller neighborhoods of Xi, if necessary), 
that the sets fli are diffeomorphic images of balls by the exponential mappings 
exp~. : B(0%.,ri) — ► fij := B(xi,ri). Consider the Riemannian manifold M :— 
Mi x • • • x M k and x :— (xi, x k ) 6 fli X • • • X O fe . Recall that if v := (vi, Vk) & 
B{Qxn r i) x ••• x B(0x k ,r k ) the exponential map exp^ is defined as exp £ (w) — 
(exp £i (vi), exp Sfc (vk)) ■ Then exp £ maps diffeomorphically the open set B(0^ ± , r\) x 
• • • x B(0 Sk ,r k ) C TM S = (TMOef, x • • • x {TM k h k onto fi x x ••• x fi fc . 

We consider the functions, defined on suitable open subsets of euclidean spaces, 

uj(t, v) := uj(t, exp s (v)), ip(t, v) := cp(t, exp £ (w)), v t ) := u t (t, exp £i (^)). 

We have that 

oj(t,vi, ...,v k ) = Ui(t,vx) H h u k (t,v k ), 

and (i, 0j) = (i, 0^, ...,0s h ) is the maximum oiui — ip. Therefore, we apply 
Theorem 8.3] to ensure, for every e > 0, the existence of real numbers bi and 
bilinear forms Bi g £j(M rai ), i = 1, k, such that 

(b l7 D Vi iP{%0s),B t ) £ V 2 ' + Ui{t,0 Si ) 
for i = 1, k, and the block diagonal matrix with entries Bi satisfies 

Bi ... 
' "A\\\ I< \ : •.. : I <A + eA z 



e 



... B k 
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where A = D 2 v ^(t, S ) and h + ■ ■ ■ + b k = ^(t, S ). Clearly 

^(t, (fe) = ^(t, x), D Vi ip(t, £ ) = D Xi <p(t, x), 

and Lemma [32] provides the equality D 2 ip(t,0s) — D 2 ip{t,x). To conclude this 
proof it remains to apply Lemma 13.71 to obtain the equivalence 

□ 



4. Comparison 

Let us state and prove our main comparison result for viscosity solutions of 
dCEEj). 

Theorem 4.1. Let M be a compact Riemannian manifold of nonnegative sectional 
curvature, and let F : Jq (M) — > R be continuous, elliptic, translation invariant and 
geometric. Let u £ USC([0,T) x M) be a subsolution and v S LSC([0,T) x M) 
be a supersolution of (|CEE|) on M. Then u < v on [0,T) x M whenever u < v on 
{0} x M. 

Proof. Since M is compact we know that M has injectivity radius %m > 0. 

Let us start noting that we may assume u and — v bounded above on [0, T) x M. 
Otherwise, for every < S < T, consider u and —v defined on the compact set 
[0, S] x M, where they are also u.s.c. and thus bounded above. Then, we apply the 
arguments of the proof to u and — v in [0, S) x M. 

Next, let us observe that for e > 0, the function u = u — ^jr-z is also a subsolution 
of u t + F(Du, D' 2 u) = on [0, T) x M. Moreover, 

(4.1) u t + F(Du,D 2 u) <-^— forDu^O, 

(4.2) u t <-j?2 for D ™ = °< and 

(4.3) lim u(t, x) = — oo uniformly on M. 

Since the assertion u < v for every e > implies u < v, it will suffice to prove the 
comparison result under the assumptions given in (4.1-4.3). 

Assume to the contrary that sup (u — v) > 0. Take / 6 T . Since M is 

[0,T)xM 

compact, u and — v are u.s.c. and (4.3) holds, we can consider for every a G N, 
m a := sup {u(s,x) - v(t, y) - af (d(x, y) 2 ) - a(t - s) 2 }, 

0<s,t<T 
x,y£M 

which is attained at some (s a ,t a , x a ,y a ) S [0,T) x [0,T) x M x M. Clearly, 

m a > sup (it — v) > 0. 

[0,T)xM 

If t a = for infinitely many a's, which we may assume are all a, then we have 
< sup (u — v) < m a = sup (u(s, x) — v(0, y) — af (d{x, j/) 2 ) — as 2 ) 

[0,T)xAf s,x,y 
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We deduce from Lemma 13.31 that lim Q ^oo af (d(x a , y a ) 2 ) = and lim Q ^oo a(t a — 
s a ) 2 — 0. By compactness, we can assume that a subsequence of (t a , s a , x a , y a ), 
which we still denote (t a , s a ,x a ,y a ), converges to a point (so, to, Xo, yo). By Lemma 
13.31 we have that xq = yo and so = £o = 0, and linia^oo m a — u(0, Xo) — v(0, Xq) = 
sup :c6M (u(0, x) — v(0, x)) < 0, which is a contradiction. 

A completely analogous argument leads us to a contradiction if s a = for 
infinitely many a's. 

Thus we may assume that there exist ao > such that s a > and t a > for 
a > a Q . By compactness and Lemma 13.31 we may also assume that x a and y a 
converge to the same point xq = yo, and in particular that x a ,y a € B{xo,r/2) for 
all a > ao, where r > is small enough such that < r < iu and conditions 
(B, C) of Section 2 hold whenever d(x,y) < r. Therefore the function d(x,y) 2 and 
hence the functions 

®a(x,y) := af(d(x,y) 2 ), tp a (s,t,x,y) := $ a (x,y) + a{t - s) 2 

are C 2 smooth on (0,T) x (0,T) x B(x ,r/2) x B(x ,r/2). 

2 — 2,+ 

Recall that V ' v(t a ,y a ) = —V ' (—v)(t a ,y a ), and if we consider the function 

*(x,y) := d(x,y) 2 

we obtain from [SJ Section 3] that 

D x ^(x a ,y a ) = -2exp 3 r o L (?/ Q ), and D y $>(x a ,y a ) = -2exp~*(x a ). 

Now we cannot directly apply Theorem 13.81 because condition (|3.ip is not 
generally satisfied due to the singularity of F (one has a serious difficulty when 

2 4- 

V ' u(s a ,x a ) contains triplets of the form (a,0,A): in this case one cannot use 
the fact that u is a subsolution to guarantee that a < C, since F((, A) — > oo as 
£ — > 0). Instead we will use Theorem 13. 1[ treating the variables s, t as if they were 
spatial variables in the stationary case, and then ignoring the information that this 
result gives about the second derivatives with respect to the variables t, s, which 
we do not need here. Bearing in mind that (s a ,t a , x a ,y a ) is the maximum of the 
function (s, t, x, y) — ► u(s, x) — v(t, y) — tp a (s, t, x, y), and setting 

A a := D (p a (s a , t a , x a , y a ), e := e a = — ■ — — , 

we obtain this way two bilinear forms P a € £ 2 (TM Xa ), and Q a e C 2 {TM ya ) such 
that 

(4.4) I ~^~ i Pa ($a, ta , Xa, ya ) : B 'x^a ($a, ta , Xa , ya ) , Pa J £ P ^{^a , Xa) , 



1)1 



(4.o) i ^ (Pa(,Sa>-> ta i^aj Va ) ; ^y^fa (^a j j ^a 7 Uoc^i Qa ) ^ ^ ^(^a 1 Ua)i 

and 

(4-6) - f- + ||A Q ||V< n )<^ + e«^. 



-Q Q 



These inequalities can be deduced from the corresponding ones in Theorem 13.11 
(just bear in mind the special form of our function cp a , and apply the inequalities 
given by Theorem l3.1l to vectors of the form (0, 0, v, w), where the zeros correspond 
to the variables s and t). 
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In our case we have 

d 

d 

• — 7^ ^Pot ('5 a ; j 7 JJa ) 2G:(t Q; S a ), 

(4.8) D x ip a (s a ,t a ,x a ,y a ) = -2af'(d(x a ,y a ) 2 )exp~^(y a ), 
-D y tp a (s a ,t a ,x a ,y a ) = 2af / (d(x a ,y a ) 2 )expy^(x a ), 

and in particular we see that 

(4.9) a a + b a = 0. 

Let us now distinguish two cases. 
Case 1. Assume that x a ^ y a - Let us consider the non-zero vectors 

Ca := -2af'(d(x a ,y a ) 2 )exp~^(y a ), 

and notice that 

L XaVa ( a = 2af'(d(x a ,y a ) 2 )expyi(x a ). 

Since u is a strict subsolution and v is a supersolution of u t + F(Du, D 2 u) = 0, 
we have that 

a a + F(( a ,P a ) <j^<0<-b a + F(L 

x a DaCa, Qa) 

(notice that here we used continuity of F off {( = 0}, and the important observation 

2.+ 

that if ((,A) e V ' u(z) and ( =/= then ((,A) is a limit of a sequence (( k ,A k ) 
with ((k,A k ) = {Dip k {z k ),D 2 ip k {z k )) for some ip k E A(F) and z k — > z). 
Thus, there is c := ^ such that 

(4.10) < c < F(L Xa ya ( a ,Q a ) - F(Ca,P a ). 
On the other hand, since F is translation invariant, we deduce 

(4.11) 0<c<F(L Xay 

a Cot i Qa) P{Ca-> Pot) — P{Cot-> -^y a x a {Qa)) P{.Coo P a ) 

Recall that A a = D 2 ip a (s a , t a , x a ,y a ) and <f> Q = af o <3>. Then 

(4.12) !>*,„*«(*> ?/) = «/'(*(*, 2/)) (£>x*(z, y), y)) = 

= -2af'(V(x, y)) (exp" 1 y, cxp" 1 x) . 

Now D 2 y ip a , the Hessian of ip a , satisfies for every vector fields X, Y on M x M 

(4.13) D^ a (M,X,y) = (V x (D<p a ),Y) = (VxK'(f)M), y) 

= a(/'($)Vx(W) + *(/'(*))£>*, y) 

- a/'($)(V x (M), y) + aI(/'(*))(M, y) 

= a/'(*)D2 iB *(A:,y) + a/"(*)X(*)<D*, F) 

= af'(<$)D 2 xy t>(X,Y) +a/"(*)(D*® L>tf)(X,y). 

In particular for every two points x,y £ M such that d{x,y) < minjiM^), *m(j/)} 
and every w <G TAfj, we consider X = y with X(x,y) = (v,L xy v) e TM X x TM y 
and we obtain 

X(*)(x,y) = Di y «f(x,y)(v,L X yv) = D x V{x, y)(v) + £>„*(&, y)(L xy v) = 0. 
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The last equality in the above expression is proved in [5] Section 3]. Therefore, if 
M has sectional curvature bounded below by some constant — Kq < 0, we obtain 
from equation (5.5) and Proposition 13.21 that 

)(v,L XaVa v) = D cp a (s ){v,L XaVa v) 

= af'(^i(x a ,y a ))D 2 ^(x a ,y a )(v, L XaVa v) 2 

(4.14) < af'(^(x a ,y a ))2K ^(x a ,y a )\\v\\ 2 , 

for every v £ TM Xa . Let us denote by Ai < • • ■ < A„ the eigenvalues of the 
restriction of A a to the subspace T> = {(v, L XaVa v) : v G TM Xa } of TM Xa x TM Va . 
The above inequality implies that Ai, A„ < 2aK , t>(x a , y a )f'(^(x a , y a ))- With 
our choice of e a , we have that 

Ai + s a X 2 < A, + — — - X 2 < X l + \X l \ < 2max{0,A„}, i = 1, ...,». 

1 +sup 1 < J <„ \Xj\ 

Since Ai + e a Af , i = 1, ...,n, are the eigenvalues of (A a + £ Q A^) \-p, this means 
that when M has nonnegative sectional curvature, that is Kq — 0, or equivalcntly 
X n < 0, we have 

Therefore, the second inequality in (|4.6[) implies P a —L Va Xa (Q a ) < {A a +e a A 2 a )\xi < 
0. Thus equation (|4.11[) . and the fact that F is elliptic imply that 

< c < F(( a , L ya Xa (Q Q )) - F(Ca, Pa) < 0, 

a contradiction. 

Case 2. If we are not in Case 1 then we may assume x a = ya for every a > uq. 
We know that 

u(s,x) - v(t,y) - af(d(x,y) 2 ) - a(t - s) 2 < u(s a ,x a ) - v(t a ,y a ) - a(t a - s a ) 2 

for all (s,t,x,y). By taking y = y a , t = t a we get that the function (s,x) i— > 
u(s,x) — af(d(x 1 y a ) 2 ) — a(t a — s) 2 has a maximum at (s a ,x a ), which (bearing in 
mind that f (0) = = /"(0)) yields 

(-2a(i a - « a ), 0,0) G P 2 '+u(s Q ,x Q ). 
Similarly, we also deduce that 

(-2a(t Q - s a ), 0, 0) g V 2 '-v{t a ,y a ). 
Since u is a strict subsolution and v is a supersolution, we get 

— £ 

—2a(t a - s a ) < — < < -2a(t a - s a ), 
a contradiction. □ 

The preceding proof can be easily modified to yield the following more general 
results. 

Remark 4.2. One can replace the compactness of M in the statement of Theorem 
14.11 by the following condition on the behavior of u and v at oo: 

(4.15) limsup u(t, x) — v(t, x) < 

(t,x) — >OG 

(this condition is meant to be empty when M is compact). 
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In the case when M does not have positive curvature, one can prove the following. 

Theorem 4.3. Let M be a complete Riemannian manifold with sectional curvature 
bounded below and positive injectivity radius. Let F satisfy conditions (A. - T)) of 
Section 2. Assume furthermore that there exist f G J~(F) and C > such that 

(4.16) tf'{t) < Cf{t) for all t > 0, 

and that F satisfies the following uniform continuity assumption with respect to the 
variable D 2 u: 

(4.17) F(C, P-5I)- F((, P) uniformly on £ P. 

Let u £ USC([0, T) x M) be a subsolution and v £ LSC{[§, T) x M) be a superso- 
lution of (|CEEp on M . Suppose that u < v on {0} x M and 

(4.18) limsupu(t, x) - v(t, x) < 0. 
Then u < v on [0, T) x M. 

Proof. Assume that the sectional curvature of M is bounded below by — Ko, with 
Kq > 0. We have, with the notation used in the proof of Theorem 14. 1[ case 1, 
following equation (|4.14[) . that A n > 0, and 

Aj + £ a ^i < 2A„ < AaK ^(x a ,y a )f'{^(x a ,y a )), 

hence 

(Ac +e a A 2 a ) (v,L Xaya v) 2 < 4aK f(x a ,y a )f'(<f(x a ,y a ))\\v\\ 2 . 
Thus, inequality (|4.6p and condition (|4.16p imply that 

P a (v) 2 - L yaXa (Q a )(v) 2 = P a (v) 2 - Q a (L Xaya v) 2 

< (A a +e a A 2 a ) (v,L Xaya v) 2 

< 4aK y(x a ,y a )f '(^(x a ,y a ))\\v\\ 2 

(4.19) <4,aK Cf(^(x a ,y a ))\\v\\ 2 . 
Let us denote 

6 a := AaK Cf{^{x a ,y a )). 

We have that hin^^oo S a = 0. From (I4.19P we obtain P a — 5 a L < L yaXa (Q a ). Then, 
equation (|4.1ip . the fact that F is elliptic, and condition (14. 17[) imply that 

< c< F{( a ,L VaXa {Q a )) - F(( a ,P a ) 

< F(( a ,P a - 5 a I) - F(C a ,P a ) 0, 

which again leaves us with a contradiction. The proof of case 2 parallels that in 
Theorem SU □ 

Remark 4.4. Condition (|4.16p is always met when one is able to take an / of the 
form 

f(t)=t k , 

with k > 2. Therefore, in the cases when F is given by the evolutions by mean 
curvature or by Gaussian curvature, (|4.16p is automatically satisfied. 
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On the other hand, condition (|4.17|) is also clearly met by the function F asso- 
ciated to the mean curvature evolution problem. Indeed, in this case the function 
A i ► F(£, A) is linear, so we have 

F((, P-SI)- F((, P) = -SF(C, I) = 6 trace - ^fj < S(n - 1), 

where n is the dimension of M. We thus recover Ilmanen's Theorem from |26j : 

Corollary 4.5 (Ilmanen). Let M be complete, with sectional curvature bounded 
below and positive injectivity radius. Let F be given by (|2 . 2[) . Let u € USC([0, T) x 
M) be a subsolution and v € LSC([0,T) x M) be a supersolution of (|MCEJ) on M. 
Suppose that u <v on {0} x M and limsup^ 2 )^ oc u{t, x) — v(t, x) < 0. Then u < v 
on [0,T) x M. 

Unfortunately, condition (|4.17p in Theorem l4.3l is not satisfied by the function F 
given by (|2.5|) corresponding to the evolution of level sets by Gaussian curvature. 
In this case, we can only apply Theorem 14. II in order to deduce a comparison result 
for manifolds of nonnegative curvature: 

Corollary 4.6. Let M be a complete Riemannian manifold of nonnegative sec- 
tional curvature and positive injectivity radius. Let F be given by ()2.5p . Let 
u G U SC( [0, T) x M) be a subsolution and v € LSC'( [0, T) x M) be a supersolution of 
(l+GCE|) on M . Suppose that u < v on {0} x M and lim. sup ^ t x ^ QO u(t,x) — v(t,x) < 
0. Then u < v on [0, T) x M . 

Given the form of the equation (|CEE|) . it immediately follows that, in all cases 
where comparison holds, one has continuous dependence of solutions with respect 
to initial data. 

Remark 4.7. Ifu,v are solutions with initial conditions g and h respectively, and 
11.9- Hl^(m) < s, then \\u - v\\ L ^ (Mx[0 ^ T)) < e. 



5. Existence by Perron's Method 

We will have to use the following estimation for the second derivative of the 
distance to a fixed point. 

Lemma 5.1. 33, p. 153] Let M be a complete Riemannian manifold whose sec- 
tional curvature K satisfies S < K < A. Suppose < r < min{iM(xo), 7r/2\/A}. 
Then, for all x £ B(xo,r) and v_LVc?(-, xq)(x), one has 

c A {d(x,x )) 2 c s {d(x,x )) , , 

— yr(v,v) < D d{;Xo)(x)[v,v) < —— rr(v,v), 

SA{d(x,x )) s s {d(x,Xo)) 
and the gradient \7d(-,xo)(x) belongs to the null space of D 2 d(-,xo)(x). 
Here s$ and c$ are defined by 

(sin(VSt))/VS, 5>0; 
s s (t) :={ t, _ 5 = 0; 

(sinh( v ^It))/v / 5, S < 0, 
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and 

cos(VSt), S > 0; 
cs(t) :={ 1, 5 = 0; 



Notice that 



cosh(V|5|t), 5<Q. 



(5.1) lim ^ = i. 

Proposition 5.2. Let F satisfy conditions (A. - T)) of Section 2, and assume 
J~(F) ^ 0. Let S be a nonempty family of subsolutions of 

(5.2) u t + F(Du,D 2 u) = 0, 
and define 

W(z) := sup{v(z) :veS}. 

Suppose that W*(z) < +oo for all z S [0, T) x M. T/jen W* is a subsolution of 
(15^1) on [0,T) x Af. 

Proof. Let G -4(-F) be such that W* — ip has a strict maximum at zq = (to,Xo). 

We may assume that W*(zo) — f{z ) = 0. 

Case 1. Suppose first that Dip(z ) ^ 0, and let us see that 

<Pt(zo) + F{D<p(z ),D 2 ip(z )) < 0. 

Define tp(t, x) := (p(t, x) + f(d(x, xoj) + (t — t ) 4 , where / G 3~(F), and observe that 

(5.3) W*(z) - i>{z) < -f(d(x, x )) - (t - t ) 4 . 

Also notice that ip £ A(F), ipt(zo) = L Pt{zo), Dip(z ) = Dip(z ), and D 2 ip(z ) = 
D 2 <p(z ). 

By definition of W* there exist z' k such that lim^oo z' k = zq and 

a k := W*(z' k ) ~ ^(4) - W*(z ) - iP(z ) = 0. 

Now, by definition of W, there exists a sequence (vk) C S such that Vk(z' k ) > 
W(z k ) — j:, which implies 

(5.4) (v k - ip)(z' k ) >a k - -. 
Since Vk < W, (|5.3[) implies 

(5.5) {v k - i>)(z) < -f{d(x, so)) - (t - t ) 4 for all z. 

Let i? be a closed ball of center Zq. Since Vk — ip is upper semicontinuous it attains 
its maximum on B at some point z k € B. From (15. 4p and (|5.5j) we get 



a fc - i < K - ^)(4) < K - VOOfe) < -f{d(x k ,x )) - (t k - t ) 4 < 0, 

and since a k — > we deduce that Zk — > zo and tk — > £o- Moreover, Vk — ip has a 
local maximum at z^. 

Since Dip(zo) = D(p(zo) ^ 0, we have Dip(z) ^ for all z in a neighborhood 
(which we may assume to be B) of Zo- Because is a subsolution and Dip(zk) =/= 0, 
we get 

V> t (z fc ) + F(Di/j(z k ), D 2 ij;{zk)) < 0. 
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Therefore, by taking limits and using the continuity of F off {C = 0} and the 
continuity oftpt, Dip, D 2 ip, we obtain 

ip t (z ) + F(D^(z ),D 2 ^(z )) = Mzo) + F(D^(z Q ),D 2 ^(z )) < 0, 

and we conclude that W* is a subsolution of (|5.2p at z . 

Case 2. Assume now that Dip(zo) = 0, and let us check that (fit(zo) < 0. Since 
if € A{F), there exist <5 > 0, lu e C(R) with w{r) = o(r), and / € ^(F) such that 

(5.6) |^(£C, - v?(z ) - Vt(z )(t - t ) | < f(d(x,x )) +uj(t - t ) 

for all z = (t,x) <E B := B(z ,5o). We may assume that cj G C 1 (R), cj(0) = = 
w'(0), and w(r) > for r > 0. Let us define 

ip(t, x) := <p t {zo)(t - t ) + 2f(d(x, x Q j) + 2w(t - to), and 

ipkfax) := <p t (zo){t - *o) + 2/(d(a;, at )) + 2w fe(* - i o), 
where (uik) is a sequence of C 2 functions on R such that cjfe — * and u)' k — > u>' 
uniformly on R. 

From (|5.6p we deduce that W* — if> has a local strict maximum at zq. On the 
other hand it is clear that (ipk ) C A(F) , and ipk — ► V> uniformly. Arguing as in Case 
1, we may find a sequence of subsolutions (vk) C 5 and a sequence of points Zk 
such that 2^ — > zo an d u fe — ^fc attains a maximum at zu- Since v k is a subsolution 
we have 

(5.7) (tpk)t(z k ) + F(Dip k (z k ) 7 D 2 ip k (z k )) < for all k, when x k ^ x , and 

(5.8) (ipk)t(zk) < 0, wn en x k = x a . 
Notice that 

(5.9) lim (ipk)t{zk) = <ft{zo)- 

k — >oo 

If Xk = xq for infinitely many fc's, we immediately deduce from (|5.8[) and (|5.9p that 
ifit(zo) < 0. 

Therefore we may assume that Xk ^ xq for all k. If we set 

Cfc :=-exp-£(x ), 
A k := D 2 d(-,x )(x k ) 
then we have that = d(x k , xq), and 

(ipk)t(zk) = </?t(-Zo) + 2w U*ft _ *o) 

= T^r/'GCfclK*, 
ICfcl 

^V fe (^) = 2/"(iai) %^ + 2/'(iai)A fc . 

Since -F is geometric we have 

(5.10) F(DMzk),D 2 Mz k )) = 2/'(|C fe |)F(^,A fe ). 

I U I 

Next, because 5 = _B(zo,<5o) is compact, we may find numbers A, 5 > such 
that the sectional curvature K of M satisfies S < K < A on B. We may of course 
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assume <5o < minjzAf (xo), 7r/2\/A}, so that we can apply Lemma [5.11 we obtain 
that Ak{(k,Ck) = 0, and for all v £ TM Xk such that wJ_Cfc we have 

CA(|Cfc|) / \ s A I \ s c *(ICk|) I \ 
777-7- (v,v) < A k (v,v) < ———(v,v). 

This implies 

^•m c A(|g|) ( T Qk®Q k \ . c s (\Ck\) ( T Ck®Qk 

13 ^(M)\ lQkT)- k -^m\ ICkT 

On the other hand, equation (|5.ip tells us that 



c^ > 1 and c J ^ < 2_ 
s A (t) ~ 2t a s (t) ~ t 
if t > is small enough. Hence we have 

sAdai) - 2|Cki w(iai) - Kfei 

for A: large enough, which we may assume are all k. By plugging these inequalities 
into (|5.11[) we obtain 



(5 ' 12) lap 

Bearing in mind that _F is elliptic and geometric, we get 

F[ ]ty Ak) - 

which combined with (|5.10l) yields 

(5.13) 2^^F(C fc ,27) < F(DMzk),D 2 Mzk)) < 2^^^, -21), 

which, thanks to condition (|2.6| . allows to conclude that 

(5.14) lim F(Dip k (z k ), D 2 ^ k {zk)) = 0. 

k — >oo 

Finally, from lf5~7)) . ([53)1 and (f5TT4|) . it follows that 

Vt(^o) < 0. 

In either case we see that W* is a subsolution of (15.21). □ 



Theorem 5.3. Let F satisfy conditions (A. - D) of Section 2, and assume that 
T(F) and comparison holds for the equation 

u t + F(Du,D 2 u) = 
u(0,x) = g(x). 

Let u andu be a subsolution and a super solution of (|5. 15[1 . respectively, satisfying 
M*(0,x) = u*(0,x) = g{x). Then w — sup{u : u < v < u, v is a subsolution} is a 
solution of (|5 . 1 5|) . 



(5.15) 
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Proof. From = < w* < w < w* < It*, we deduce that w*(0, x) = w(0,x) = 
w*(0,x) = g(x). On the other hand w* is a subsolution by Proposition 15. 2[ and 
w* < u by comparison, hence w* = w by definition of We claim that w* is 
a supersolution. This implies w* < by comparison, hence = w = w* and 
consequently w is a solution. 

Let us prove the claim. Suppose to the contrary that w* is not a supersolution. 
Then there exist zq — (to,Xo) and a C 2 function tp such that (w* — p)( z ) > — 
(w* — <p){zo) for all z, and either 

(5.16) (pt(zo) + F(Dtp(z ),D 2 (p(z )) < 0, when Dtp(z ) ^ 0, or 

(5.17) <Pt(zo) < 0, when Dtp(z ) = 0. 

By replacing tp(t, x) with the function tp(t, x) + d(x, Xq) 4 + (t — to) 4 on a neigh- 
borhood of zq we can furthermore assume that 

(5.18) (ic - <p)(t, x) > d(a?, x ) 4 + (i - t ) 4 - 
Let us denote 

C/ 5 := {(t, z) : dfo x ) 4 + (t - i ) 4 < £ 4 }- 
Case 1. In the case when (|5.16[) holds, by continuity of tp t , D<^, D 2 ^> and F, we 
can find r > such that 

<^ t (z) + F(^(z),^ 2 <^(z)) < 

for all z £ U2r, that is tp is a subsolution on Ui T , and obviously the same is true of 
<p := tp + r 4 /2. 

From (|5. 18[) we have that 

4 ^4 

(5.19) w(z) > w*( z ) - y > </?(2) + y for all z £U 2r \ U r . 
Now let us define 

W(z) = { max {<P( z )> w ( z )}' if z e U r\ 

1 otherwise. 

By using Proposition 15.21 and equation (|5.19|) , it is immediately checked that W is 
a subsolution. We have W = w outside B(z , r), but 

(5.20) swp(W - w) > 

because, by definition of w*, there exits a sequence x n )} converging to (to, xq) 
such that lvaxw(t n ,x n ) — w*(to, xo), and consequently we have 

Hm(W(t n , x n ) - w(t n ,x n )) > lim(^(t„, x n ) - w(t n , x n )) = r 4 /2 > 0. 

On the other hand, W(0,x) = w(0,x) — g(x), because we could of course have 
taken r > small enough such that (0, x) £ B(zq, r). We deduce (from comparison 
again) that W < u and consequently W < w, which contradicts (|5.20[) . 

Case 2. On the other hand, in the case when (|5.17p holds, since tp £ A(F), there 
exist 5 > 0, u £ C(K) with w(r) = o(r), and / G T(F) such that 

\tp(x,t) - tp(z ,to) - tp t (t ,x Q )(t - t )\ < f(d(x,x )) +u(t - t ) 

for all z = (t,x) £ B := B(z ,8o). We may assume that w £ C 1 (R), lj(0) — = 
uj'(0), and w{r) > for r > 0. Let us define 

ip(t, x) = (p(z ) + tp t (z )(t - t ) - 2f(d(x, x )) - 2u>(t - t ). 
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Then w^ — tf) attains a strict minimum at zq. Also notice that Dip(z) ^ for z =/= zq. 
Arguing as in Case 2 of the proof of Proposition [521 one can show that 

(5.21) lim F{D^{z),D 2 ^(z)) = 0. 

z^z 

By combining this with the continuity of tpt and the fact that tpt(zo) = Vt( z o) < 0, 
we can find an r > such that 

Mz) + F(D^(z),D 2 ^(z)) < 

for all z 6 f/ 2r , z ^ zq. The rest of the proof is identical to that of Case 1 (just 
replace ip with ?/>). □ 



Let us now show how to apply the above Theorem in order to construct solutions 
of (|5.15p . We will need to use the following stability result. 

Lemma 5.4. Assume that u.s.c. (respectively l.s.c.) functions Uk are subsolutions 
(super solutions, respectively) of (|CEE|) . Assume also that {uk} converges locally 
uniformly to a function u. Then u is subsolution (super solution, respectively) of 
(ICEEl) . 

Proof. Suppose that ip £ A(F) and u—ip attains a strict local maximum at (to, Xq). 
The convergence of the subsolutions Uk allows us to find a sequence of local maxima 
(tk, %k) of Uk — (p which converges to (to, xq). Then, by a similar argument to that 
of the proof of Proposition 15. 2\ one can show that u is a subsolution of (|CEEj) at 
(t ,x )- □ 



Let M be a compact Riemannian manifold. Assume that comparison holds for 
the equation (|5.15p . Let us first produce solutions of (|5.15|) for initial data g in the 
class A(F). 

Let us define 

u(t, x) = —Kt + g(x), and u(t, x) = Kt + g(x), 

where K := sup^g^- \F(Dg(x), D 2 g(x))\ (which is finite because g 6 A(F ) and M 
is compact). It is immediately seen that u is a subsolution and u is a supersolution 
of (|5.15p . and obviously u*(0,x) = u*(0,x) = g(x). According to Theorem 5.3 and 
comparison, there exists a unique solution u of (|5.15|) . 

Now take g a continuous function on M. According to Proposition 12. 8[ we can 
find a sequence gk of functions in A(F ) such that gk — ► g uniformly on M. Let Uk 
be the unique solution of (|5.15|) with initial datum gk- By Remark 14. 7\ (uk) is a 
Cauchy sequence in C([0, oo) x M), hence it converges to some u € C([0, oo) x M) 
uniformly on [0, oo) x M. Then by Lemma 6.1 it follows that u is a solution with 
initial datum u(0,x) = g(x). 

Therefore we can combine this argument with Theorems 14.11 and 14.31 to obtain 
the following corollaries. 

Corollary 5.5. Let M be a compact Riemannian manifold of nonnegative sectional 
curvature, g : M -tla continuous function, and let F : Jq (M) — > R be continuous, 
elliptic, translation invariant and geometric. Then there exists a unique solution of 
(fCEE) on [0,oo) x M. 
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Corollary 5.6. Let M be a compact Riemannian manifold, g : M — ► R a continu- 
ous function, and let F satisfy conditions (A.-H) of Section 2. Assume furthermore 
that (|4.16p and (|4.17p are satisfied. Then there exists a unique solution of (|CEE|) 
on [0, oo) x M. 

Corollary 5.7 (Ilmanen). Let M be a compact Riemannian manifold, g : M — » 
R continuous. There exists a unique solution u of the mean curvature evolution 
equation (jMCEJ) such that u(0,x) = g(x). 

Corollary 5.8. Let M be a compact Riemannian manifold, g : M — ► M continu- 
ous. There exists a unique solution u of the positive Gaussian curvature evolution 
equation (|+GCEj) such that u(0,x) = g(x). 

Corollary 5.9. Let M be a compact Riemannian manifold, g : M — ► R continu- 
ous. There exists a unique solution u of the mean curvature evolution equation in 
arbitrary codimension (given in Example \2.3\) such that u(0,x) = g{x). 

When M is not compact, analogous corollaries can be established if one addi- 
tionally demands that the initial datum g be a (positive) constant outside some 
bounded set of M, and that i(M) > 0. The proof is similar (replacing uniform 
convergence on M with uniform convergence on compact subsets of M) . 



6. Geometric consistency & Level set method 

Theorem 6.1. Let 9 : R — > K be a continuous function, and let u be a bounded 
continuous solution of (|CEE[) . Then v = 9 o u is also a solution. Moreover, if 9 
is nondecreasing and u is a subsolution (resp. super solution) then v — 9 o u is a 
subsolution (resp. super solution) as well. 

Proof. Assume first that 9 is monotone. We may consider a sequence of smooth 
functions 9k with nonvanishing derivatives, converging uniformly to 9 over the 
bounded range of u. Hence by Lemma 15.41 we may directly assume that 9' ^ 0. 
Notice that g = 9~ x satisfies g' ^ too. 

Suppose first that 9' > 0. Let tp E A(F) and assume that 9ou — (p attains a local 
maximum at zq. If we denote ijj = g o ip, it is not difficult to check that ip € A(F), 
and u — ip clearly attains a local maximum at (to,xo). Consequently 

i> t {z )+F{Dip{z ),D 2 ip{z Q )) <0 

if Dtp(to,Xo) ^ 0, and ^t(toj^o) < otherwise. But Dip(zo) — if and only if 
Dip(zo) = 0, and 

Di> = g'(ip)D<p 

D 2 ip = g"{p)Dip ®D%p + g'{ip)D 2 ip. 
Since F is geometric and <?' > 0, one immediately sees that 

<pt(zo) + F(Dp>{z ),D 2 p>{z )) = , ( ) (zq) (^Oo) + F{Di,(z ), D 2 yj(z ))) < 

if Dip(zo) 0, and (pt(to,xo) = g i/Ji Z0 \\ V'tC^o) < otherwise. This shows that 9ou 
is a subsolution. 
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If 9' < 0, the same argument tells us that if u is subsolution (respectively super- 
solution), then v is supersolution (respectively subsolution). In order to establish 
the result for continuous functions, it is enough to observe that a continuous func- 
tion can be uniformly approximated by a sequence of locally monotone functions. 
Then a local application of Lemma [5T4l yields the result. □ 

Now one can show that, if comparison and existence hold for (|CEEp (e.g. when 
M is a compact Riemannian manifold of nonnegative curvature), then for every 
compact level set To there is a unique, well-defined, level set evolution Tt of To by 
the geometric curvature evolution equation corresponding to (jCEEj) . 

Let g be a continuous function on M with ro = {x G M : g{x) — 0}, and assume 
that Tq is compact. We may also assume that g is constant outside a bounded 
neighborhood of To, and in particular bounded. Let u be the unique solution of 
flCEE} with u(0, •) = g. We define 

r t = {x G M : u(t,x) = 0}. 

Theorem 6.2. Assume that comparison and existence hold for (|CEE|) . Let g : M — » 

R be a continuous function satisfying To — {x G M : g(x) = 0} and such that g 
is constant outside a bounded neighborhood of To. Let u be the unique continuous 
solution of (|CEEj) with initial condition g. Then 

r t = {i£ M : u(t,x) = 0}. 

Proof. This is a consequence of Theorem 16.11 and the comparison principle. It 
follows exactly as in the case M — R n , see [13l Theorem 5.1], or [17l Chapter 4], 
for instance. □ 

Corollary 6.3. The definition ofT t = {x G M : u{t, x) = 0} does not depend upon 
the particular choice of the function g satisfying Tq = {x G M : g(x) = 0}. 

It can also be checked that the evolution To i— » IC(t)To := T t thus defined has 
the semigroup property 



Some other properties of the evolutions can be established as in the case M = R n . 
For instance, in the case of the evolution by mean curvature, it is possible to show 
that if Tq = dU is a smooth connected hypersurface with positive mean curvature 
with respect to the inner unit normal field, then T t continues to have positive mean 
curvature as long as it exists, in the sense that 



(which admits a unique viscosity solution, see [5]). 

However, one has to be very cautious and not take it for granted that all the 
usual geometrical properties of the generalized evolutions by mean curvature or 
by Gaussian curvature could be immediately extended from the Euclidean to the 
Riemannian setting. As a matter of fact, many of these properties are very likely 



IC(t + s)=JC(t)IC(s). 



T t = {xeM : v(x) = <}, 
where v is the solution of the stationary problem 




on U; 

on T = dU, 
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to fail in the case of manifolds of negative curvature. We will present several 
counterexamples and related conjectures in Section 8. 

7. Consistency with the classical motion 

In this section we suppose that equation ([CEEp arises from a classical geometric 
evolution for hypersurfaces in M. We establish the consistency of the level set 
evolution equation with this classical geometric motion. 

More precisely, suppose (r*) tg [oT] ^ s a f amir y of smooth, compact, orientable 
hypersurfaces in M evolving according to a classical geometric motion, locally de- 
pending only on its normal vector fields and second fundamental forms. In partic- 
ular, we shall assume that T t is the boundary of a bounded open set Ut C M and 
that there exists a family of diffeomorphisms of manifolds with boundary 

0* : Fo -» U~ t , t £ [0, T] , 

such that: 

(i) 0° =Id, and, 

(ii) for every x € To the following holds: 

(7.1) (x) =G(y (t, </>* (*)) , V T v (t, </? (x))) , 

where v (t, •) is a unit normal vector field to Tt, and the linear map 

V r v (t, x) : {TT t ) x 3^Vj V (t, x) G (TT t ) x 

and V T stands for the orthogonal projection onto (TTt) of covariant de- 
rivative in M. 

Classical motion by mean curvature corresponds to taking / (u, V r f) = tr (— V r ^) 
whereas classical motion by Gaussian curvature is defined by / (u, V r ^) = det ( — V v 
The level set evolution equation induced by (|7.1[) is of the form (|CEE|) where F is 
related to G through formula (|2.1[) . As before, we assume that F is elliptic, trans- 
lation invariant and geometric. In this case we already know that F is continuous, 
and in fact smooth off {£ = 0}, because F is of the form (|2.1j) with G smooth. 

Define d : [0, T] x M — ► R, the signed distance function from T t , as: 

r dist(x,r t ) if xe u t 
1 ' x) :_ \ -dist 0, r t ) if x e m \ u t . 

Lemma 7.1. There exist constants K, 6 > such that d is smooth in 
Is := {(t,x) e [0,T] x M : \d(t,x)\ < 6} 

and 

(7.2) \d t + F(Dd,D 2 d) \ <K\d\ in I s . 
Proof. We consider geodesic normal coordinates from T t , 

$ (t, s, x) := exp x (sv (t, x)) , 

assuming that v(t,x) points towards the interior Ut for every x € T t . Clearly, for 
s small enough, 

d(t, $ (t, s,x)) = s. 

Given xq S Tq there exists a neighborhood V of xq in Tq and an interval (— r, r) 
such that $ (t, •, •) is a diffeomorphism from (— r, r) x $ (V) onto its image X t := 
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$((-?-, r) x 0* (V)) for t £ [0,T]. Note that $ is also smooth in t. Denote by 
\& (t, •) the inverse of $ (i, •, •) and write 

*(t,y) := (p(i,y),X(t,y)). 

Now for x £ (f) 1 (V) and s € (— r, r) we have X (t, $ (t, s, ar)) = a; and p (t, $ (t, s, x)) — 
s. Both X and p are smooth in t, and clearly p = d in IJtero Tl {0 x Xt = h- 

In order to prove (|7.2|) it suffices to note that, since, for x € T 4 , Dd(t,x) = 
u(t,x)^0, 

r(t,x) :=d t +F(Dd,D 2 d) 

is a smooth function vanishing for x £ T t . This gives (|7.2|) locally; a global bound 
then follows by the compactness of IV □ 

Next we state and prove the main result of this section. 

Theorem 7.2. Let u be the unique viscosity solution to the level set equation hCEBfl 
with initial datum u\t=a — d\t=o- Then, for every t £ [0, T], the zero level set of 
u(t,-) coincides with T t : 

T t = {x £ M : u (t, x) = 0} . 

Proof. Define 

v (t, x) := e tK ((d (t, x) V 0) A 5/2) , 
where K is the constant given by (|7.2p . We shall prove that v is a viscosity su- 
persolution to equation (|CEE|) . Clearly, v\t=o > u\ t= o A 5/2, and, by Theorem 16. 1[ 
it A (5/2) is a viscosity solution to (|CEE[) as well. The comparison principle (Theorem 
14. lj) then will ensure that v > u A {8/2). In particular, 

(7.3) {x £ M : u (t, x) > 0} C U t , t £ [0, T] . 

On the other hand, we shall prove that 

w (t, x) := e~ tK ((d (t, x) A 0) V (-5/2)) 

is a viscosity subsolution to (|CEEp . Noww| t= o < u\t=o V(— 5/2), and the comparison 
principle will imply that w < uV (—5/2). This, together with (|7.3[) yields 

{a; eM : u (t, x) > 0} = ?7 t , t £ [0, T] . 

Now take e > and let u e := u + e; this is again a viscosity solution to (|CEE() . It 
turns out that 

v E (t, x) := e tK (((d (t, cc) + e) V 0) A 5/2) , 
w e (t, x) := e-* K (((d (t, x) + e) A 0) V (-5/2)) , 
are respectively super and subsolutions to (|CEE|) provided e is much smaller than 
5 - namely, small enough to ensure that v £ and w £ are smooth in the regions 
< d(t, x) + e < 5/2 and —5/2 < d(t,x) +e < 0, respectively. Applying the 
comparison principle as we did before we ensure that 

{x £ M : u (t, x) > -e} = {x £ M : d (t, x) > -e} . 

Letting e go to zero we conclude that the zero level set of u (t, ■) is precisely T t , as 
we wanted to prove. We now show our claim that v is a supersolution to (|CEEI) . 

Start noticing that Lemma [7 . 1 1 ensures that, for t £ [0,T] the following holds: 
(7.4) 

v t (t, x)+F (Dv (t, x) , D 2 v (t, x)) > 0, if < d (t, x) < 5/2, 
v t (t, x) > 0, if d (t, x) < or d (t, x) > 5/2. 
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(7.7) 



Let (to, xq) G [0, T] x M and ip G A (F) be such that v — ip has a local minimum at 
(to,Xo). Without loss of generality, we may assume that 

(7.5) ip(t Q ,x Q ) =v(t ,x ), 
and that 

(7.6) if < v, locally around (to,xo). 

Since the level sets d(t,x) = c e (—6,5) are smooth hypersurfaces, necessarily 
d (to, xo) 7^ 5/2. If moreover (to,xo) T t , then (|7.5|) . (|T.(3|) and the smoothness of 
v imply that, at (to, xo) one has tp t = v t and Dip = Dv. Therefore, using (|7.4[) we 
conclude that: 

ip t (to, xo) + F (Dp (t , x ) , D 2 ip (t , x )) > 0., if Dip (t , x ) ^ 0, 
tpt (to,xo) > 0, otherwise. 

Now we shall prove that the above identity also holds when (to,xo) G IV Let 
Q := Ute[o t] W x ^V this ' s precisely the set of zeroes of d. As \Dd\ = 1 on Ig, 
we infer that Q is a smooth hypersurface of [0, T] x M; since (|7.5[) and (|7.6p imply 
that v — ip has a minimum at (io,xo), the following holds for every tangent vector 
(t,£) G TQ (totXQ y. 

(7.8) ^(*o s xo)r + D^(td,i )(f) = 0; 

moreover, for every curve 7 : (—1, 1) — > L t[) with 7 (0) = xq, and 7' (0) = £ we have, 

(7-9) ^(io,7(*))| s =o<0. 
From identity (|7.8[) we deduce that 

(7.10) (tp t (t ,x ) , Dp (t , xo)) = A (d t (to, x ) , Dd (to, x )) , 
for some A G K, whereas (|7.9j) merely states that: 

(D 2 p (to, xo) C, < - (V^ (to, x ) , 7" (0)) . 
Taking into account that (Vd (to, 7 (0)) , 7' (0)) = 0, we obtain: 

(7.11) (D 2 tp (to, xo) £, £} < - (Vy> (i , xo) , 7" (0)) 

= -A (Vd (t , xo) , 7" (0)) = A (D 2 d (t , xo) £, . 

Given a smooth curve 77 : (—1,1) — > Q such that 77 (0) = (io,x ) and rf (0) = 
— (dt (to, Xo) , Vd (to, xo)) the following holds: 



d 
dt 



ip(ri (t)) \ t=0 = -A (dt (t ,x ) 2 + \Dd(t ,xo)\ 2 ^ 



Since necessarily v(rj(t)) = v (rj (0)) = for t > sufficiently small, (|7.5j) and (|7.6p 
imply that A > 0. 

Now, (H21) trivially holds if A = 0. Suppose A > 0, then using ([TTU]) . f7~TT]) we 
deduce 

-F (Dp (t , x ) , D 2 ip (t Q , x Q )) < -F (XDd (t , ar ) , XD 2 d (t , x )) 

= -XF (Dd(to,xo),D 2 d(to,x Q )) , 

and 

Vt (to, x ) = Ad t (t , x ) = -XF (Dd (t ,x ) , D 2 d (to, x )) . 
Therefore, we conclude that <p satisfies ()7.7|) at (to,xo). 
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A completely analogous proof shows that v £ is a supersolution, and w and w £ 
are subsolutions to (|CEE|) . □ 

Remark 7.3. In the very special case of the evolution by mean curvature in arbi- 
trary codimension (see Example 1 2. 3 P the above proof breaks down because the sets 
T t are no longer hypersurfaces of M, but fc-codimensional submanifolds. The con- 
sistency of the generalized motion (given, e.g., by Corollarv l5.9|) with the classical 
evolution thus remains an open problem (whose solution would probably require a 
careful analysis of the properties of the distance functions to the submanifolds T t , 
and of the eigenvalues of their Hessians, similar to the study carried out in [2] for 
the case M — W l ). 



8. Counterexamples and conjectures 

In this final section we provide some counterexamples showing that many well 
known properties of the evolutions by mean curvature in W 1 fail when M is a 
Riemannian manifold of negative curvature. 

Example 8.1. When (M,g) is the Euclidean space equipped with the canonical 
metric, Ambrosio and Soner have proved in [2] that the distance function \d\ is 
always a supersolution to the mean curvature equation (|MCE|) . This is no longer the 
case of a general Riemannian manifold, as the following example shows. 

Let (M,g) be a surface of revolution embedded in ]R 3 , locally parameterized by: 

/ x (-7r,7r) 3 (s,6) i ► (r (s) cos 9, r (s) sin 0, s) G K 3 , 

where I C M is an open interval and r S C°° (I) with r > p > 0. In these 
coordinates, the metric takes the form: 

r> (s) 2 + 1 
r(sf 

Suppose G I and r' (0) = 0; then the "Equator" s — is a geodesic of M, denote 
it by Tq. The classical evolution by mean curvature starting from To is constant 
in time. Therefore, the corresponding signed distances satisfy d(t, •) = d(0, •) = d 
for every t g R (we shall assume d > for s > 0). Let us next explicitly compute 
d. The geodesies of M that are orthogonal to Tq are of the form (s (t) , 9) with 
9 G (—it, tt) constant. Take such a geodesic and assume that is parameterized by 
arc length and satisfies s' > 0. In particular, since its tangent vector (s',0) must 
be of norm one, 

(8.1) v(s(t)) 2 (s'(t) f = l, v(s) := y/r' (,s) 2 + 1. 

Clearly, d (s (t) , 0) = t and d s d (s (t) , 9) s' (t) = 1. Identity (|8.1[) and our assump- 
tion s' > allow us to conclude that d s d (s, 9) — v (s) and Dd = (1/v, 0). Finally, 
a direct computation gives: 

/ Dd \ _ r' 
\\Dd\) ~ rv' 

The function \d\ will be a supersolution to the mean curvature equation provided 
r' (s) sign (s) < 0. This is always the case if the curvature of M remains nonncg- 
ative. On the other hand, \d\ will a subsolution whenever the curvature of M is 



\Dd\ div 
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nonpositive everywhere. Finally, taking for instance r (s) = 1 + cos 2 s we are able 
to produce a \d\ that is not a supersolution, neither a subsolution to the mean 
curvature equation. 

Conjecture: If M has nonnegative sectional curvature then \d\ is always a super- 
solution. If M has negative curvature then there always exists To such that \d\ is 
not a supersolution. 

Example 8.2. When M = R", Evans and Spruck p31 Theorem 7.3] showed that if 
To, Tq are compact level sets and T t , T t are the corresponding generalized evolutions 
by mean curvature, then 

dist(r ,f ) < dist(r f ,f t ) 

for all t > 0. This result fails for manifolds of negative curvature. For instance, let 
M = {(x,y,z) el 3 : x 2 + y 2 = 1 + z 2 } be a hyperboloid of revolution embedded 
inM 3 . Let 

r = {(*, y,z) eM:z = 0}, 

and 

t = {(x,y,z)GM:z=l}. 

Then 

r t = T for all t > 0, and dist(r , f ) > 0, 

but 

dist(r t , f t ) = dist(r , f t ) -> as t -> oo. 

Conjecture: Evans-Spruck's [T3l Theorem 7.3] result holds true for all manifolds 
of nonnegative sectional curvature, but fails for all manifolds of negative curvature. 

Example 8.3. In the case M = M™ it is well known that equation (jCEEp preserves 
Lipschitz properties of the initial data. Namely, if g is L-Lipschitz and u is the 
unique solution of (|CEE|) then u(t, •) is L-Lipschitz too, for all t > 0; see [TT1 
Chapter 3]. Since the proof of Theorem 7.3 in |13j remains valid for any manifold 
provided that one assumes the Lipschitz preserving property of (|CEEp , the preceding 
example also shows that (|CEE|) does not preserve Lipschitz constants when M is a 
hyperboloid of revolution. 

Conjecture: The equation (jCEEj) has the Lipschitz preserving property if and only 
if M has nonnegative sectional curvature. 



9. Appendix: Existence and uniqueness of viscosity solutions to a 
(nonsingular) general parabolic equation 

In this appendix we present the standard definition of viscosity solution and state 
existence and comparison result for viscosity solutions to non-singular parabolic 
fully nonlinear equations. 
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Definition 9.1. Let M be a finite-dimensional Riemannian manifold, and a con- 
tinuous function F : (0, T) X M X M X J 2 Af — > R. Consider the parabolic equation 

(9.1) u t + F(t, x, u, Du, D 2 u) = 0, 

where u is a function of (t, x). We say that an USC function u : (0,T)xM-tl is 
a viscosity subsolution of the partial differential evolution equation provided that 
a + F(t, x, u(t, x), C, A) < for all (t, x) £ (0, T) x M and (a, C, A) £ V 2 ^u{t, x). 
Similarly, a viscosity supersolution of (|9.1| is a LSC function it : (0, T) x M — > R 
such that a + F(t, x, u(t, x), C, A) > for every (t, x) € (0, T) x M and (a, C, A) £ 
V 2 '~u(t, x). If u is both a viscosity subsolution and a viscosity supersolution of 
u t + F(t, x, u, Du, D 2 u) = 0, we say that u is a viscosity solution. 

Remark 9.2. If u is a subsolution of u t + F(t, x, u, Du, D 2 u) = and F is lower 

2.+ 

semicontinuous, then a + F(t, x, u(t, x), C, A) < for every (a,£,A) <E V u(t,x) 
and every (t, x) £ (0,T) x M. A similar observation applies to supersolutions when 
F is upper semicontinuous, and to solutions when F is continuous. 

Theorem 9.3. Let M be a compact and F : (0,T) x M x K x J 2 Af -> K &e 
continuous, proper, and such that 

(1) there exists 7 > with 

7 (r - s) < F(t, 2;, r, C, Q) - F(t, x, s, (, Q) 

for all r>s, (t, x, (, Q) € (0, T) x M x J 2 M); 

(2) there exists a function to : [0, 00] — > [0, 00] with lim t ^ + ^(t) — and such 
that, for every a > 0, 

F(i,y,r,ae3cp~ :1 (a;) 1 Q) - F(t, x, r, -aexp^T 1 ^), P) < w (ad(a;,y) 2 

/or all t e (0, T), x, y £ M, r£R,P£ £ 2 (TM X ),Q £ £ 2 {TM y ) with 

-(^ + ||4j)l< ( _° Q )<Ax + ea^ 

and e a = (2 + 2||v4 ct ||)~ 1 7 where A a is the second derivative of the func- 
tion (p a (x,y) = ^d(x,y) 2 at the point (x,y) £ M x M with d(x,y) < 
mm{i M {x), iM{y)}- 

Assume that u £ USC([0,T) x M) is a subsolution and v £ LSC([0,T) x M) is a 
supersolution of (|9.ip on M and u < v on {0} x M. Then u <v on [0, T) x M. 

The proof of this result is a combination of the ideas in the proof of Theorem 
8.2] with the new techniques for second order nonsmooth analysis on manifolds 
developed when dealing with the stationary case in [3 Theorem 4.2]. See also 
Theorem 13.81 which should be used in this proof. We leave the details to the 
reader's care. 

Analogous results to Corollary 4.6 and Corollary 4.8 stated in [5] can also be 
obtained in a similar way for these general evolution equations. We say in such 
cases that "comparison holds" . 

As it is customary in these cases, whenever the Comparison Theorem holds there 
is no difficulty in applying Perron's method (see [HI HE]) to show that the problem 

/ u t + F(t,x,u(t,x),Du(t,x),D 2 u(t,x)) = 0, 
1 j I u(0,x)=g(x). 
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has a unique bounded viscosity solution u on M, provided that comparison holds 
and one is able to find a viscosity subsolution u and a viscosity supersolution u 
such that u t (0, x) = u*(0, x) = g{x). In fact u is given by 

W(t, x) = sup{u>(t, x) : u< w <u and w is a subsolution of (|9.1[) }. 
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